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We study neutrino spin-flavor oscillations in presence of rapidly varying external fields. The 
general formalism for the description of neutrino oscillations in arbitrary rapidly varying external 
fields is elaborated. We obtain the effective Hamiltonian which determines the evolution of the 
averaged neutrino wave function. We apply the general technique to neutrino oscillations in rapidly 
varying magnetic fields. The special case of the constant transversal and twisting magnetic fields is 
studied. We evaluate the effect of neutrino spin-flavor oscillations in rapidly varying magnetic fields 
of the Sun. The numerical solutions of the Schrodinger equation with the Hamiltonian accounting 
for the constant transversal and twisting magnetic fields are presented. We compare them with our 
approximate analytical transition probability formula and reveal good agreement at high frequencies 
of the twisting magnetic field. 
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I. INTRODUCTION 

For the first time the idea of neutrino oscillations was 
theoretically predicted in Ref. 0. Since then many ex- 
perimental and theoretical studies of neutrino oscilla- 
tions have been carried out. One of the most interest- 
ing problems in neutrino physics is the solar neutrino 
deficit. Nowadays it is experimentally established (see, 
e.g., Ref. 0) that the disappearance of solar electron neu- 
trinos can be accounted for by the LMA-MSW solution 
0, 0| • The process of neutrino oscillations is likely to be 
the most reliable explanation of the solar neutrino prob- 
lem. There are several theoretical models of solar neu- 
trino oscillations such as mentioned above LMA-MSW 
solution. However other scenarios like spin-flavor preces- 
sion (see Refs. [H, IS 13 ) are also considered. 

Recently resonant neutrino spin-flavor oscillations were 
studied in Ref. [Tol ] where an attempt was made to repro- 
duce the data of the major solar neutrino experiments 
using the peak profiles of the solar magnetic field. Dif- 
ferent approach to the solar neutrino problem was made 
in Ref. jllj. Supposing that neutrino spin-flavor preces- 
sion played a subdominant role in comparison with flavor 
oscillations the combined action of these two mechanisms 
was examined. It was also possible to impose the con- 
straints on the characteristics of a neutrino, namely to re- 
strict its magnetic moment, and the strength of the solar 
magnetic field. Solar electron neutrino transitions into a 
non-electron antineutrino were studied in Ref. [T^ . On 
the basis of the resent SNO data the indication of the neu- 
trino Majorana magnetic moment was obtained in that 
paper. A global analysis of the solar neutrino problem 
solution with help of the spin-flavor precession was given 
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in Ref. |13|. Two flavor oscillations scenario in the opti- 
mized self-consistent magneto-hydrodynamical magnetic 
field profile was adopted. A careful analysis of the pre- 
dicted solar neutrino fluxes was performed in Ref. Q 
using up to date helioseismological data. 

We studied neutrino oscillations in electromagnetic 
fields of various configurations in our previous works. 
First it is necessary to mention that the Lorentz invari- 
ant formalism for the description of neutrino spin-flavor 
oscillations was elaborated in Ref. It allows one to 
study neutrino oscillations in arbitrary electromagnetic 
fields. For example, neutrino oscillations in an electro- 
magnetic wave and in a longitudinal, with respect to the 
neutrino momentum, magnetic field were considered in 
Ref. ^5|- I n R e f- we worked out the method for the 
investigation of the neutrino evolution equation solution 
near the resonance point. This technique is of great im- 
portance when the evolution equation cannot be solved 
analytically. In Ref. ^3] we have applied this method to 
the studying of the parametric resonance in neutrino os- 
cillations in amplitude modulated electromagnetic wave. 
The major feature of Refs. 0,EJ was the use of the per- 
turbation theory in solving the neutrino evolution equa- 
tion. It should be noted that the perturbation theory was 
also used in Ref. |l8j| where the parametric resonance in 
neutrino oscillations in matter with periodically varying 
density was studied. In Ref. we proposed the Lorentz 
invariant quasi-classical approach for the description of 
neutrino spin oscillations in arbitrary non-derivative ex- 
ternal fields. This method was applied in Ref. [2(j to the 
description of neutrino spin oscillations in gravitational 
fields. 

In this paper we study neutrino oscillations in pres- 
ence of general rapidly varying fields. It should be noted 
that the influence of rapidly varying fields on mechanical 
oscillations was investigated in Ref. |2l|. It is also well 
known that there are certain analogies between mechani- 
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cal and neutrino oscillations (see, for instance, Ref. |22|). 
In Sec. m we start from the neutrino evolution equation 
with the Hamiltonian which accounts for the neutrino 
interaction with rapidly varying fields. Note that we do 
not fix the explicit form of the Hamiltonian. Then we 
derive the new effective Hamiltonian governing the time 
evolution of the averaged neutrino wave function. There- 
fore the obtained new Hamiltonian allows one to study 
the neutrino conversion in presence of arbitrary rapidly 
varying external fields. Our result is also beyond the 
perturbation theory because no assumptions about the 
smallness of rapidly varying fields are made. Then, in 
Sec. IIIII we apply the general technique to neutrino os- 
cillations in rapidly varying magnetic fields. We investi- 
gate the combination of constant transversal and twisting 
magnetic fields. The neutrino conversion in solar mag- 
netic fields is discussed. In Sec. II VI we present the numer- 
ical solutions of the Schrodinger equation for the neutrino 
system interacting with constant transversal and twisting 
magnetic fields and compare them with the approximate 
analytical solutions found in the present work. 



II. GENERAL FORMALISM 

Let us consider the evolution of the two neutrinos v = 
(yi, V2), which can belong to different flavors and helicity 
states. The evolution of the system is described by the 
Schrodinger type differential equation, 



■ du tt 



(1) 



where the Hamiltonian H involves external fields, e.g., 
electromagnetic field, interaction with matter and the 
characteristics of the neutrinos. Here we do not spec- 
ify the explicit form of the Hamiltonian but just suppose 
that it is decomposed into two terms, 



H = H +H, H(t + T) = Hit). 



(2) 



The first term in Eq. Q - i?o, corresponds to the neu- 
trino interaction with constant or slowly varying exter- 
nal fields. In presence of only this term the solution of 
Eq. can be easily found. The solution is known to 
be periodical with the typical frequency £!o ~ 1/L e f{, 
where L e g is the oscillations length. The second term in 
Eq. - Ti(t), corresponds to rapidly varying external 
fields. The frequency u> — 2tt/T should be much greater 
than f2o : w 3> f2o- For example, we can suppose that 
the components of the Hamiltonian depend on time like 
Tt{t) ~ cosujt or smut. Note that we do not make any 
assumptions about the strength of varying external fields. 

We will seek the solution of the Eqs. QJ and J5J in the 
form (see also Ref. plf 1 ). 



(3) 



In Eq. iJSJ the function is the small rapidly oscillating 
one with zero mean value. The mean value / of a func- 
tion f(t) is the time averaging over the period T. The 



function vo(t) in Eq. (0 is the slowly varying one during 
T. Substituting Eq. © in Eq. JJJ we obtain 



(4) 



There are two groups of the terms in Eq. : rapidly and 
slowly varying ones. In order for this expression to be an 
identity the terms of each group in the left-handed side 
of the equation should be equal to the terms of the same 
group in the right-handed side. 

First let us discuss rapidly varying terms. The func- 
tion £ standing in the left-handed side belongs to this 
group. Despite the function £ is small, its derivative can 
be, in principle, not small. Indeed, it is proportional to u> 
and thus £ is the great quantity. There are three rapidly 
varying terms in the right-handed side of Eq. Q, namely 
Hvo, T~t£, and H ^. However the second and the third 
terms in this group are the small ones because they con- 
tain the small factor £. Thus we can drop the terms H£ 
and Hq^ when we discuss rapidly varying terms. Finally 
one receives the equation for the function £ 



(5) 



It should be noted that the terms in Eq. containing 
£ can contribute to the phase of the wave function if 
long evolution times are considered. However, when we 
discuss possible applications of this approach in Sec. IIIII 
only one period of the transition probability variation 
is taken into account. This problem is also analyzed in 
Sec. |EI Equation © can be easily integrated and we 
find its solution 



H(t)dt v (t). 



(6) 



Here we take into account that variations of the func- 
tion vq are small. Therefore Eq. © presents the time 
dependence of the function £. 

Now let us proceed to the description of Vq evolution. 
We can omit i£ and Ttvo in Eq. |J2J| because we have 
already proved that these terms are equal [see Eq. (JSJ]. 
Averaging residual terms in Eq. over the period T 
and allowing for 



i> = !>o, H i/ = H v , H £ = H £, = 0, 
we obtain the equation for the function vq 

dl/ ° it , nTc 
i-r- = H u + H£. 

at 



(7) 



We cannot neglect the term 7i£ in Eq. Q since it is the 
product of two periodically varying functions. The mean 
value of this quantity is not equal to zero. With help 
of Eq. © equation can transformed into the more 
common form 



■ dp ° IT 
at 



(8) 
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where 



Constant transversal magnetic field Bo, 



H eS = H -iH (J H&t). 



(9) 



In derivation of Eqs. iJSJ and 10 we do not make any 
assumptions about the smallness of the interaction de- 
scribed by the Hamiltonian H. Thus our result is beyond 
the perturbation theory used in Refs. [l6L IrH Il8| . 

The Hamiltonian H can be always represented in the 
form, 

H=(vh), 

where h is the arbitrary 3D-vector. Using this represen- 
tation we discuss two cases: 



• Twisting magnetic field B(r). 

Note that if a magnetic field is constant in time we may 
consider only its transversal component with respect to 
the neutrino velocity. These magnetic field configura- 
tions were studied separately in the majority of works 
devoted to neutrino spin-flavor oscillations. For instance, 
only constant transversal component of a magnetic field 
was taken into account in Refs. [UEl, and the effect of 
only twisting magnetic field on neutrino oscillations was 
considered in Refs. 0,S0. 

In our case the Hamiltonians Hq and H have the form 



• H o = (erj_ ■ hi), H = a 3 h 3 , 

• H = (<r± ■ hi), H = a 3 h 3 . 

Here we introduce 2_D-vectors cr± = (01, 02) an d hi = 
(hithz). 

In the former case the application of the developed 
in the present paper technique leads to no interesting 
consequences. Indeed, using Eq. JJ]J we find that 



ffeff = Oi • hi) - ilh 3 (jh 3 dt), 

where 1 is the unit matrix. It is well known that terms in 
the Hamiltonian proportional to the unit matrix can be 
omitted. For example, if we study neutrino oscillations 
in matter with rapidly varying density n(t), it will cause 
no essential effect on the oscillations process. However, if 
one accounts for all components of the vector h in H, i.e. 
H = (<r • h), it could result in some interesting effects. 
This case requires special careful examination. 

Now we consider the latter case. First we introduce 
two auxiliary 2D- vectors a = hi and b = J adt. Again 
using Eq. © we obtain that 



_ (V/2 vBo \ 



and 



H = 







/iBe 




(10) 



(11) 



where /i is the neutrino magnetic moment. In Eq. 110|) 
we introduced the quantity, 



V 
~2 



Am 2 
4E 



-e 



Gf 
V2 



n c ff, 



vac 
<— > 



where O is the function of the vacuum mixing angle 6 
(the explicit form of for various transitions of the vn, 
VjR type can be found in Ref. |!j), Am 2 is the difference 
of the neutrino mass squared, E is the neutrino energy 
n e ff is the effective matter density, to is the frequency of 
the transversal magnetic field variation, Gf is the Fermi 
constant. 

Using Eqs. (P|)- (|llfl it is possible to derive the expres- 
sion for -Hetf, 



H cS =0-3/13 - it (a. ■ b) + 0-3(01^2 - a 2 6i) 



03 h 3 + (ai&2 - CLzbi) 



Therefore the off-diagonal rapidly varying terms can shift 
the resonance point. It was shown in Ref. ,19j that both 
axial- vector and tensor and pseudotensor interactions can 
cause the neutrino spin precession. Thus rapidly varying 
external fields of the mentioned above types, e.g., elec- 
tromagnetic fields, interaction with moving and polarized 
matter etc., will nontrivially affect the neutrino oscilla- 
tion process. We study below one the possible examples. 



H, 



efit 



V/2 - {{jlBY/lo LiB 

fiB -V/2 + {hB) 2 /uj 



(12) 



It should be noted that the parameters B and Bq must 
have non zero values. However it is interesting to discuss 
the case Bq — > 0. When a neutrino interacts only with 
twisting magnetic field the problem of neutrino oscilla- 
tions can be solved exactly (see, e.g., Ref. 0)- One can 
explicitly write the formula for the transition probability 



III. NEUTRINO SPIN-FLAVOR OSCILLATIONS 
IN RAPIDLY VARYING MAGNETIC FIELDS 

In this section let us discuss the neutrino evolution in 
matter under the influence of a combination of the two 
types of magnetic fields 



where fi' = y/(V - w) 2 /4 + (fiB) 2 . It can be shown that 
in the limit of rapidly varying twisting magnetic field the 
transition probability approaches to zero. Indeed 



uj^>fiB,V 



-t 



Thus we can see that the transition probability vanishes 
when a neutrino interacts only with rapidly varying twist- 
ing magnetic field. This result also follows from Eq. (jl2|l . 

Since Bo and B do not depend on time in Eq. (|12Jl , we 
can solve Eq. JSJ for the case of the effective Hamiltonian 
given in Eq. 1|12[) . The transition probability is expressed 
in the following way, 



Pit) = Asm. — 



Tit 



L 



where 



A = 



and 



[V/2 - {pB)*/uY + (fiBo) 2 



l = V= y/[V/2 fcB)7w] a + (a^o) 2 



(13) 



(14) 



(15) 



Thus the validity of the developed method (the fre- 
quency of the twisting magnetic field u> should be much 
greater than the characteristic "frequency" of the sys- 
tem at the absence of this field) can be expressed as the 
constraint, 



uj » il - ^/(V/2) 2 + (pB ) 2 . 

Now let us discuss the resonance conditions in neu- 
trino oscillations in rapidly varying magnetic fields. From 
Eqs. I|13|) and i|14|l it follows that if the condition is sat- 
isfied, 



V _ (nBf 



(16) 



then A ~ 1 and the transition probability can achieve 
great values. This phenomenon is analogous to resonance 
amplification of spin-flavor oscillations. However, in our 
case the resonance is attained not due to the zero value 
of the parameter V. 

We again emphasize that the strength of the twisting 
magnetic field B can be not small in the proposed tech- 
nique. The only assumption made consists in the great 
value of u>. Therefore we can discuss two cases: 

• B » B, 

• B > B a . 

In the former case the value of A in Eq. I|14|l is great inde- 
pendently of B. This case corresponds to the usual spin- 
flavor precession in constant transversal magnetic field. 
However, in the latter case the value of Ao = A(B = 0) 
is much less than unity. Hence the transition probability 
at the absence of the additional twisting magnetic field is 
small. In this case we can choose the parameters so [see 
Eq. I |lti|l ] that the amplitude of the transition probability 
is great. It is this situation which is of interest because 
the essential neutrino spin-flavor conversion is attained 




FIG. 1: Neutrino flux propagation in the magnetic fields of 
the Sun, Bo and B are the vectors of constant tranversal and 
twisting magnetic fields. 



under the influence of compact perturbations of a mag- 
netic field against a rather weak constant background 
magnetic field. Such magnetic field configurations can be 
created experimentally or observed in some astrophysical 
media. Magnetic field in perturbations can have twist- 
ing structure. From Eqs. (|14f) - i|ltty one can derive the 
restriction imposed on the B and Bo, 



B 
~B) 



» 1. 



(17) 



Now we consider the possible application of the de- 
veloped technique to spin-flavor neutrino oscillations in 
magnetic fields of the Sun. The solar magnetic field is un- 
likely to be only either constant transversal or twisting. 
Therefore we can apply the method elaborated in this 
paper to the description of the solar neutrino conversion. 
The propagation of the neutrino flux in the combina- 
tion of constant transversal and twisting solar magnetic 
fields is schematically depicted in Fig.^ We consider one 
of the possible channels of neutrino oscillations, namely 
v e h conversion. First we should estimate the pa- 

rameter V in Eq. (|10|) . In our case the function O and the 
effective matter density have the form (see, e.g., Ref. 0), 



e = 



1 + cos 26> v 



1 



Let us discuss a neutrino with the following properties: 
Am 2 ~ 10" 5 eV 2 , 6> vac - tt/4, which are not excluded 
by the modern experimental results (see, for instance, 
Ref. |23). We take the neutrino energy E ~ lOMeV, 
which corresponds to the 8 B solar neutrinos. Matter is 
supposed to consist mainly of the hydrogen, i.e. n n ~ 0, 
n e — n p (as a consequence of the system electroneutral- 
ity) and has the density d ~ 1.4 g • cm~ 3 , which is close 
to the mean density of the solar matter. For these pa- 
rameters we obtain that V/2 ~ 10~ 15 eV. 

Now let us evaluate the strength of the magnetic fields, 
Bq and B, necessary for the 10% conversion of the ini- 
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tial i/ e L beam. We suppose that neutrinos conversion 
occurs along the distance D = t ~ 9.2 x 1O~ 2 -R ~ 
3.2 x 10 14 eV _1 , where Rq is the solar radius. We also 
suppose that neutrinos have the transitional magnetic 
moment [i = 10~ n /iB, where /ib is the Bohr magneton. 
Setting P(t) = 0.1 in Eq. (Q3J) we obtain B ~ 17.6 kG. 
From the condition (|17|l we can derive the strength of the 
twisting magnetic field: B ~ 3.2i?o — 56.3 kG. The res- 
onance condition i|16|l is satisfied if u ~ 1.0 x 10 _14 eV. 
The close values of the twisting magnetic field frequencies 
(3 x 10~ 15 eV) were considered in Ref. Q- 

Background matter perturbations can cause twisting 
magnetic fields generation. If a neutrino interacts with 
solar plasma, the frequencies of these magnetic fields can 
be evaluated with help of the magnetic hydrodynamics 
(MHD) methods. We suppose that the twisting magnetic 
field is induced by a MHD wave with the wave vector per- 
pendicular to Bo- The frequencies of such MHD waves 
can achieve values up to 1.5 x 10 5 eV for B ~ 17.6 kG 
(see, for instance, Ref. |24()- The frequency of the twist- 
ing magnetic field, used in the present work, satisfies this 
condition. It should be noticed that solar plasma fluctu- 
ations were studied in Ref. 251 . 



IV. NUMERICAL SIMULATIONS 

The analysis carried out in the previous section shows 
that the influence of rapidly varying magnetic fields on 
neutrino oscillations is very important in various astro- 
physical environments, for example, in v e L v^r oscil- 
lations in magnetic fields of the Sun. In order to sub- 
stantiate the correctness of the approach developed in 
this paper for the description of this problem we obtain 
the numerical solutions of Eq. Q with the Hamiltonian 
given in Eqs. (|10l) and {TTJ. The numerical solutions 
account for all terms in Eq. Q. Then we get the transi- 
tion probabilities on the basis of the numerical solutions. 
We also compare these numerical transition probabilities 
with our approximate formula for the transition proba- 
bility [Eq. i|13fl ] which is derived in the limit of rapidly 
varying Hamiltonian H.. The numerical simulation is car- 
ried out with MATLAB language for technical computing 
(version 6.5). 

The measurable quantity of neutrino oscillations is the 
transition probability. The probabilities of the neutrino 
oscillations are shown on Figs. 12141 versus r = 1.5£. We 
also compare them with the approximate analytical for- 
mula (JT3J). We remind that we take V/2 = 10~ 15 eV, 
Bo = 17.6 kG, B = 56.3 kG. The resonance frequency 
is uj = 1.0 x 10 _14 eV. The magenta line corresponds 
to the transition probability computed with help of the 
numerical solution of Eq. for the taken parameters. 
The black line represents the transition probability given 
inEq. 

It is noticed in Sec. [H] that the small factors 7i£ and 
i?o? can contribute to the phase of the neutrino oscil- 
lations. However if we consider only one period of the 
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FIG. 2: Neutrino transition probability for v e L <-> v^b. oscil- 
lations, lu — 1.0 x 1CP 14 eV. The magenta line corresponds 
to the transition probability found with help of the numerical 
solution of Eq. JQ, the black line represents the transition 
probability given in Eq. jH) . 
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FIG. 3: Neutrino transition probability for u £ l *-* v^b. oscil- 
lations, ui = 2.0 x 10 -14 eV. The magenta line corresponds 
to the transition probability found with help of the numerical 
solution of Eq. l(T). the black line represents the transition 
probability given in Eq. 1131 . 



neutrino transition probability variation, then there is 
rather good agreement between numerical and approx- 
imate analytical solutions. Moreover it clearly follows 
from Figs. 1 2141 that the phase shift decreases if ui be- 
comes greater. Therefore it can be seen in Figs. 12141 that 
the numerical expression for the transition probability 
approaches to the approximate formula given in Eq. i|13|) 
at great frequencies of the twisting magnetic field. This 
comparison proves the validity of the elaborated tech- 
nique. 

At the end of this section we note that the transition 
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FIG. 4: Neutrino transition probability for u e L *-* v^b, oscil- 
lations, u = 4.0 x 10~ 14 eV. The magenta line corresponds 
to the transition probability found with help of the numerical 
solution of Eq. Q, the black line represents the transition 
probability given in Eq. I|13p . 

probability presented in Fig. |5] corresponds to the pa- 
rameters that almost coincide with the resonance ones 
(see Sec. 1111(1 . Despite of the rather low frequency (u> = 
1.0 x 10~ 13 eV), there is a good agreement between the 
predicted oscillations parameters such as amplitude of 
the transition probability [Eq. 11411 ] as well as oscillations 
length [Eq. I|15|l ] and those obtained from the numeri- 
cal solution. The best parameters coincidence happens 
within the first period of transition probability variation. 

V. CONCLUSION 



of the averaged neutrino wave function. This new ef- 
fective Hamiltonian enabled one to study neutrino os- 
cillations in arbitrary rapidly varying external fields. It 
is worth mentioning that, in contrast to Refs. [Til Il7|. 
the strength of rapidly varying fields has not been lim- 
ited. Therefore the elaborated method was beyond the 
perturbation theory since we have not carried out any ex- 
pansions over the strength of the external fields. We have 
examined the resonance conditions in the case of general 
rapidly varying external fields. It has been revealed that 
matter with rapidly varying density caused no effect to 
neutrino oscillations. We have demonstrated that rapidly 
varying axial- vector (e.g., interaction with moving or po- 
larized matter), tensor and pseudotensor external fields 
could result in nontrivial effects in neutrino oscillations. 
On the basis of the derived new effective Hamiltonian we 
have considered neutrino spin-flavor oscillations in the 
combination of constant transversal and twisting mag- 
netic fields. The new effective Hamiltonian for this sys- 
tem has been obtained. We have analyzed various limit- 
ing cases (Bo ^> B and B ^> Bq). The limit Bq — > has 
also been discussed. It has been shown that described 
mechanism of neutrino oscillations could be important 
for the neutrino spin-flavor conversion in magnetic fields 
of the Sun. We have evaluated the strengths of magnetic 
fields and frequency of the twisting magnetic field nec- 
essary for the 10% neutrino conversion from z/ e j, to v^r. 
The numerical solutions of the exact (without neglecting 
any terms) Schrodinger equation for the two neutrino 
system interacting with constant transversal and twist- 
ing magnetic fields have been found. We have compared 
them with the approximate analytical formula for the 
transition probability. The very good agreement at high 
oscillations frequencies has been established. 



In conclusion we note that in this paper we have 
examined neutrino oscillations in general rapidly vary- 
ing fields. We have started with the usual Schrodinger 
type evolution equation. The Hamiltonian involved both 
slowly and rapidly varying in time terms. Instead of solv- 
ing the evolution equation directly we have derived the 
new effective Hamiltonian which described the evolution 
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